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Abstract 

Quantum inequalities are constraints on how negative the weighted average of the renormalized 
stress-energy tensor of a quantum field can be. A null-projected quantum inequality can be used 
to prove the averaged null energy condition (ANEC), which would then rule out exotic phenomena 
such as wormholes and time machines. In this work we derive such an inequality for a massless 
minimally coupled scalar field, working to first order of the Riemann tensor and its derivatives. We 
then use this inequality to prove ANEC on achronal geodesics in a curved background that obeys 
the null convergence condition. 
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I. INTRODUCTION 


In general relativity it is possible to have exotic spacetimes that allow snperlnminal travel, 
closed timelike curves, or wormholes, so long as the appropriate stress-energy tensor is 
available. Although general relativity does not provide any restrictions on quantum field 
theory does. These constraints are called energy conditions or quantum (energy) inequalities. 
The simplest energy conditions are bounds on projections of the stress energy tensor at each 
point in spacetime, but those are easily violated by quantum fields, even by free fields in fiat 
spacetime. But by averaging, we can produce conditions that are not so easily violated. A 
quantum inequality bounds an average of T^y over a localized part of a timelike path, and 
an averaged energy condition bounds the energy along an entire geodesic. 

A good technique to rule out exotic phenomena jH is to use the achronal averaged null 
energy condition (achronal ANEC), which requires that the null-projected stress-energy 
tensor cannot be negative when averaged along any complete achronal null geodesic. 


Tabtt > 0 


( 1 ) 


where 7 is an achronal null geodesic (also called a null line), i.e., no two points of 7 can be 
connected by a timelike path, and £“ is the tangent vector to 7. 

Ref. proved ANEC for null geodesics traveling in fiat spacetime (though there could 
be curvature elsewhere) using a quantum inequality. In previous work | 3 | we studied ANEC 
in the case of a classical curved background, meaning a spacetime generated by matter that 
obeys the null energy condition. 


Tabtt>Q, 


( 2 ) 


at all points and for all null vectors We conjectured a particular form for a curved-space 
quantum inequality, and from that we were able to show that that a quantum scalar field 
in a classical curved background would obey achronal ANEC. Here we complete the proof 
by demonstrating a curved-space quantum inequality (somewhat different from the one we 
conjectured before) and using it to prove the same conclusion. 

The rest of the paper is structured as follows. In Sec. [TTl we state our assumptions and 
present the ANEC theorem we will prove. We begin the proof by constructing a parallel¬ 
ogram which can be understood as a congruence of null geodesic segments or of timelike 
paths, as in Ref. [^. In Sec. Hill we present and discuss the general quantum inequality 
of Fewster and Smith . Secs. HVR^ apply that general inequality to the specific case 
needed here, using results from our previous application of Fewster and Smith’s inequality 
in Ref. . In Sec. IVIII we present the proof of the ANEC theorem of Sec. [TTl using the 
quantum inequality. Finally, Sec. IVIIII is a summary of our results and discussion of some 
open problems. 

We use the sign convention (—,—,—) in the classification of Misner, Thorne and Wheeler 
[^. Latin indices (in small or capital letters) from the beginning of the alphabet will denote 
all coordinates; those from the middle of the alphabet will denote only spatial coordinates. 
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II. THE THEOREM 


A. Assumptions 


We consider a spacetime M containing a null geodesic 7 with tangent vector and define 
a “tubular neighborhood” M' around 7 , which is composed of a congruence of null geodesics 
as in Ref. [^. 

Then we define Fermi-like coordinates on M' as follows j^. First pick some point p on 
the geodesic 7 . Let E(^u) = and pick a null vector at p such that = 1, and two 

unit spacelike vectors E(^x) and E(^y) at p, perpendicular to E(^u) and E(^^) and to each other, 
giving a pseudo-orthonormal tetrad. Then the point q = (m, n, x, y) in these coordinates 
is found by traveling unit distance along the geodesic generated by vE(^^) + xE(^x) + yE{y), 
parallel transporting , and then unit distance along the geodesic generated by uE(^u) ■ 

We suppose that the curvature inside the tubular neighborhood M' obeys the null con¬ 
vergence condition, RabV°‘V^ > 0 for any null vector V. This will be true if the matter 
generating this curvature obeys the null energy condition, Eq. (|2]). 

We require that in M' the curvature is smooth and obeys the bounds. 


\R 


abed 


< Rr. 


( 3 ) 


and 


\R, 


abed,a 


< RL 


\R. 


abed,a/3 


<r: 


R 


abed,a/3'y 


<r 

^ -^max 5 


( 4 ) 


in the coordinate system described above, where the greek indices a,/3,'y,... take values 
v,x,y but not u, and Rmax,-Rmax)-^max)-^max finite numbers but not necessarily small. 
These bounds need not apply outside M'. 

Finally, we consider a quantum scalar field in M. Inside M' it is masseless, free, and 
minimally coupled but outside M' we allow interactions and different curvature couplings. 
For further details see Sec. II E of Ref. [^. 


B. The theorem 

Theorem 1. Let {M,g) be a spacetime and 7 an achronal null geodesic and suppose that 
around 7 there is a tubular neighborhood M'. We suppose that the curvature is bounded in 
the sense of Sec. Ill Al and the causal structure of M' is not affected by conditions outside M' 
jsl . Let Tab be the renormalized expectation value of the stress-energy tensor of a minimally 
coupled quantum field in some Hadamard state co. 

Then the ANEC integral, 

/ OO 

dXTabti\TiX)) (5) 

-OO 

cannot converge uniformly to negative values on all geodesics F(A) in M'. 


C. The parallelogram 

We will use the {u,v,x,y) coordinates of the Fermi-like coordinate system defined in 
Sec. IllAl Let r be a positive number small enough such that whenever 1^1, |a;|, \y\ < r, the 
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point (0,n,a;, I/) is inside the tubular neighborhood M' dehned in Sec. Ill A[ Then the point 
{u,v,x,y) G M' for any u. Dehne the points 


<h(M,n) = (M,n, 0,0), 

with V hxed and u varying. Write the ANEC integral 

POO 

A{v)= / duTuu{^{u,v)). 


( 6 ) 


( 7 ) 


As in Ref. we suppose that, contrary to Theorem 1, Eq. o converges uniformly to 
negative values, and show that this leads to a contadiction. 

Given any positive number vq < r we can End a negative number —A greater than all 
A(v) with V G (—VqjVo). By uniform continuity, it is then possible to find some number Ui 
large enough that 


ru+iv) 


' U— ( r ) 


for any v G (—Vq, Vq) as long as 


duTuu{^{u,v)) < -A/2, 


■u+(n) > Ml 
■u_(m) < —Ml 


( 8 ) 


(9a) 

(9b) 


We define a sequence of parallelograms in the (m, v) plane, and integrate over each par¬ 
allelogram in null and timelike directions. The parallelograms have the form 


V G (-Mo, Mo) 

M G (m_(m),m+(m)) , 

where m_(m),m+(m) are linear functions of v defined below.. 

Let / be a smooth sampling function supported only on (—1,1) and normalized 


(10a) 

(10b) 


daf{aY = 1 • 


( 11 ) 


'-1 


Then we can take a weighted integral over the whole parallelogram. 


rvQ 


dv fiv/vof 


f'U+iv) 


’ -VQ 


' U— ( r ) 


duTuui^(u,v)) < -VoA/2. 


( 12 ) 


We choose a velocity V and define the Doppler shift parameter 

V 1 - 1 / 

We pick any hxed number a with 0 < a < 1/3 and let 

to = d-^r, 

and choose 

Mo = to/(\/2(5). 


(13) 

(14) 

(15) 
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FIG. 1. The parallelogram <!>(«, u), v G (— u G {u-{v),u+{v)), or equivalently 
t G {-to, to), T] G (-??o,%) 


Then as 1/ —)■ 0, 5 —)• cxo and to, Vq —>■ 0. We define 

r/o = ui + to^/V2 (16a) 

u±(v) = ±r]o + ■ (16b) 

The points 

<S>v(v,t) = ^(v+^,^') , ( 17 ) 

with |? 7 | < ijo and |t| < to, are the same parallelogram described above, but parameterized 
in a different way (see Fig. [T]). For constant rj the paths are timelike and in ffat space 
parametrized W proper time. In curved spacetime t is approximately the proper time as 
shown in Ref. [3|. 

Now we change variables in Eq. flT^ using the Jacobian 


d{u,v) 1 
d{ri,t) 


(18) 


to get 



dtTuu{^v{'n,'t))f{i/iof < 


-Ato/2 . 


(19) 


We will show that this upper bound conflicts with a lower bound that we will derive using 
quantum inequalities on the paths given by fixing r] and varying t in ^v{vA)- 
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III. A GENERAL QUANTUM INEQUALITY 


Quantum inequalities are bounds on weighted averages along a timelike path of projec¬ 
tions of the stress-energy tensor Tab- The general form is 



dtf(t)T,t{w{t))V‘‘V^>-B, 


( 20 ) 


where w{t) is a timelike path parametrized by t, V is a. vector field onto which the stress- 
energy tensor will be projected, f{t) is a smooth sampling function, and B is some positive 
number depending on the choice of quantum field, the spacetime, the projection direction V, 
and the function /. In this paper we will apply the general quantum inequality of Fewster 
and Smith to the case of Tuui^v) appearing in Eq. (IT^ . 

Following Refs. [3, [^ , we define the renormalized stress-energy tensor, 

x^x' 


The quantities appearing in Eq. fl^ are dehned as follows. The operator is the point- 
split energy density operator, 

Taf = Va ® - gab'g^'^'Vc O Vrfi , (22) 


which is applied to the difference between the two point function and the Hadamard series. 


H{x,x') = 




-F + ^Vj{x,x')a^^{x,x')\n +^Wj{x,x')a^{x,x') 

(^+[x,x) V f 


(23) 

We have introduced a length I so that the argument of the logarithm in Eq. fl23|) is dimen¬ 
sionless. The possibility of changing this scale creates an ambiguity in the definition of H, 
but this ambiguity for curved spacetime can be absorbed into the ambiguity involving local 
curvature terms discussed below j3]. For simplicity of notation, we will work in units where 
I = 1 . 

In the first term is the Van Vleck-Morette determinant, and a is the squared invariant 
length of the geodesic between x and x', negative for timelike distance. In flat space. 



a{x, x') = —r]ab{x — x')°‘{x — x'Y ■ 

(24) 

By F((T+), for some 

function F, we mean the distributional limit 



F(a+) = lim F(a,), 

£^■0+ 

(25) 

where 


ae{x, x') = (t(x, x') + 2ie{t{x) — t{x')) + . 

(26) 


In some parts of the calculation it is possible to assume that the two points have the same 
spatial coordinates, so we define 

T = t-t', (27) 

and write 

F(a^) = F(T.)^\imF(T.), (28) 
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where 


(29) 


Tf = T — le . 


The Hadamard series can be written 


H{x,x') = ^ Hj{x,x '), 
i=-i 


(30) 


where the snbscript j shows the power of a in the term. Following the notation of Ref. j^, 
we let denote the snm of all terms from H_i throngh Hj. 

The qnantity Q is added “by hand” to ensnre that the stress-energy tensor is conserved 
Bnt since we will be interested here in projection on a nnll vector i, Q will not contribnte, 
becanse = 0. 

The term Cab handles the possibility of inclnding local curvatnre terms with arbitrary 
coefficients in the definition of the stress-energy tensor, 
terms inclnde 


From Ref. 10 we find that these 


^^'>Hab = 2R.ab - 2gabaR - gabRV2 + 2RRab (31a) 

^'^^Hab = R-ab — '^Rab “ gab^R/2 — gabR’^‘^Red/2 + 2R’^'^Racbd ■ (31b) 

So we must include a term in Eq. fl?I|) given by a linear combination of Eqs. f|31ap and 
flHbl) . However, we keep only first order in R, and ignore those terms that vanish on null 
projection, so for our purposes. 

Cab = a ^^^Hab + b ^^^Hab ~ 2aR^ab - b{R^ab - '^Rab) , (32) 

where a and b are undetermined constants. 

From Ref. we have the definition 

H{x,x') = ^[H{x,x') + H{x',x) +iE{x,x')] , (33) 

where iE is the antisymmetric part of the two-point function. We will let Ej be the part of 
E involving a\ define a “remainder term”, 

j 

Rj = E — ^ ^ Ek , (34) 

k=-l 

and let 

Hj{x, x') = - [Hj{x, x') + Hj{x', x) + iEj{x, x')] (35a) 

H(j) {x, x') = ^ [R(j) (a;, x') + Rq) {x\ x) + iE{x, x')] . (35b) 

We will use the Fourier transform convention 

pOO 

f{k) or f^[k] = / dxf{x)R^^ . (36) 


7 








We can now state the qnantnm ineqnality of Ref. j^, on a timelike path w{t) with the 
stress-energy tensor contracted with nnll vector held 

n VT 


dTg{ty{riXr)Uw{t)) > 


{g^gWTX'tfHX] (-^,0 


-t- / dtg\t)Cabtt, 


(37) 


where g{t) is a smooth fnnction with compact snpport and the operator 6* denotes the 
pnllback of the fnnction to the path, 




The snbscript (5) means that we inclnde only terms throngh j = 5 in the snms of Eq. ([23]). 
However, as we proved in Ref. j^, terms of order j > 1 make no contribntion to Eq. fimi . 
Thns we can write Eq. fl57|) with w(t) = f) ^ specihc valne of g and the stress energy 

tensor nnll contracted with vector held pointing only in the u direction with 7“ = 1, 

/ CO 

dr(?(f)2(Tr)(Mt))>-R, (39) 

■OO 


with 

B= ^F(-e,0-/ dtgyt){2a + b)X^, (40) 

J 0 ^ J —OO 

where 

= g{t)g{t')TS^HX'^{t),w{t')) , (41) 

F denotes the Fonrier transform in both argnments according to Eq. fl36l) . and we nsed the 
fact that Ruu = 0 according to Ref. j^. 


IV. CALCULATION OF 


We will now evalnate Eq. fITUD in the case of interest. In this section we will calcnlate 
and thus F{t,t'). In Sec. jVl we will Fourier transform F{t,t'), and in Sec. IVTI we 
will hnd the form of B in terms of limits on the curvature and its derivatives. 

To simplify the calculation we will evaluate in a coordinate system [t, x, y, z) 

where the timelike path w{t) points only in the t direction, the z direction is perpendicular 
to it, and x and y are the previously dehned ones. More specihcally t and z are 


5 ^u + 5v 5 — 5v 


(42) 


where we extend the dehnition of t from Sec. HI] to cover the whole spacetime. The new null 
coordinates u and v are dehned by 


u = 


t + z 

'TT 


and are connected with u and n. 




V = 6v. 


(43) 


(44) 
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The operator can be written 



( 45 ) 


If we define C, = z — z' and u as the u coordinate of x, the center point between x and x', 
we have 



A. Derivatives of 


For the derivatives of H^i it is simpler to use Eq. (05]). We have 



(47) 


In flat spacetime it is straightforward to apply the derivatives to H_i. However in curved 
spacetime, there will be corrections first order in the Riemann tensor to both a and its 
derivatives. 

We are considering a path w whose tangent vector is constant in the coordinate system 
described in Sec. Ill Al The length of this path can be written 



where Ax = x — x' and x” = x' + AAx since dx°‘/dX is a constant. 

Now a is the negative squared length of the geodesic connecting x' to x. This geodesic 
might be slightly different from the path w. However, this deviation results from the con¬ 


nection, which is hrst order in the curvature (times the coordinate distance from the origin 
— see Eq. (9) of Ref. [^). Thus the distance between the two paths is hrst order, and the 
difference in the metric is second order in the curvature (see Eqs. (25,27) of Ref. 0 )- The 
difference in length in the same metric due to the different path between the same two points 
is also second order. All these effects can be neglected, and so we take a = — s^. 


Now using Ref. we can write the hrst-order correction to the metric. 


(49) 


9ab Vab “ 1 “ “ 1 “ -^ba 7 

where Fab is given by Eq. (29) of Ref. because the hrst step foTx = y = 


0 is in the v 


direction and the second in the u direction. By the symmetries of the Riemann tensor the 
only non-zero component is 



(50) 


where we took into account the diherent sign conventions. Putting this in Eq. fl48l) gives 



( 51 ) 
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So to first order in the curvature, 


a{x, x') = —s{x, x'Y = —— 2 / (iXFifyAx^Ax"". 

Jo 

We define the zeroth order a, 
and the first order, 

(j^^\x^x') = —2 [ dXFfjyAx'^Ax'^ 

Jo 

= -2/ dX [ dK{l-K)Ryuyu{Kx''^,x''^)x''^x''^Ax^Ax^ 


'0 JO 
rl ri 


=-2 / dX dy{(i - y)Ryuvu{,y,x ^)Ax'^Ax'", 

Jo Jo 

where we dehned ^ = x”'" and changed variables to y = ni. Now to hrst order. 


a 


( 1 ) 


(T+ ’ 


and the derivatives. 


d d \ f 1\ 


dx^ dx'"^ J \a+ J 


= ± + W(.) _ 

A ' A ^ 


C=o 


rl r® 


2^2 / 


( 1 ) ^( 1 ) 


U ,U' 


( 1 ) 


- — to-v.. 


r 


4 ,uu 


Now we can take the derivatives of a , 


=-2 I dXX^ I dy{i - y)Ryuyu{y,x"")Ax"Ax 


d 


=-2 dXX / dyRuyuy{y,x”")Ax"Ax" . 


Similarly, 


a« = -2/ dX{l-X)^l dy{l-y)Ryyy^{y,x"")Ax^Ax^ 


d 


=-2 dX{l-X) / dyRuyuy{y,x''")Ax"Ax". 


For the two derivatives of 

= -2 t x"")Ax^Ax^ . 

Jo 

Now we can assume purely temporal separation, so Ax" — Ax" — Tj\Pl and 


x"= (i"+ =,("_=), 

\/2 


(52) 

(53) 

(54) 


(55) 


(56) 


(57) 


(58) 

(59) 

(60) 
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( 61 ) 


where z = (z + 2:')/2 and t" = t' + Ar. Then the derivatives of are 

£-2 / /•! 

rpSplit Tj _ 

>0 Jo 

- 2 V 2 [ dX{l-2X) [ dyRuyuy{y,x'"’)T 


1 = 4 ^2^4 ( 4-12 / dX dy{i-y)Ryuyu{y,x"'’) 


>0 


>0 


+ / dXil-X)XRyy^yix”^,x"^)R 
Jo 


Let ns dehne the locations = {kx'^,x^) and 


< = + z)R'’ - z). 


Then Eq. fl6Tl) can be written 


rv-iSplit Tj- 

-l-uu' -^-1 - 


1-2 




dX 


12 / 

Jo 


+2\/2(l-2A) / - (1 - A)Ai?si)^iii(x")r^ 


The derivatives of iL_i can thns be written 

tS‘H-1 = V" 


r \ TT 


T ( A + 2/l(^i'^) ) + + ■\y3{t,T) 


T_ 


T_ 


where the |/j’s are smooth fnnctions of the cnrvatnre, 

yi{t,T) = f dXYi{t") y 2 (t,T)= [ dX{l-2X)Y2(t") yzii.r) = f dX{l- 
Jo Jo Jo 

with 


(62) 


(63) 


(64) 


■A)Ayi(«"), 

(65) 


r.(i") = 


n((") = 


27r^ 

1 r 


2'k‘^ 


dK{l- Y){t'' + zfRyyyu{x''^) , 

) 

dhi(t T z'^ Ryyyyi^X , 


1 


Ys{t") = -—R^,^,{x "), 

where x" and x" are dehned in terms of t" by Eqs. fl60|l and 


(66a) 

(66b) 

(66c) 


B. Derivatives with respect to r and u 

Ref. calcnlated .ff(i), bnt for points separated only in time. Let us use coordinates 
(T, Z, X,Y) to denote a coordinate system where the coordinates of x and x' differ only in 
T. Ref. a gives 

^(i)(T, T') = H_,{T, T') + Ho{T, T') + H^{T, T') + Uri{T, T') , (67) 


^(o)(T, T') = T') + H,{T, T') + -iR,{T, T '), 


( 68 ) 
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where 


^_i(T,T') = 

Ho{T,r) = 


47r2(T-T'-ie)2 ’ 


487r2 


Rtt{x) - -R{x) In (-(T - T' - zef) 


^i(T,r) = 


(T-T 

6407r2 


•i \2 r 


-Rtt,tt{x) — -UR[x) 


-5 ( ORim + ifl, tt(S) ) In (-(T - r - ^,f) 


The order-0 remainder term is 


Ro{T,T') = 


327r2 


dn{ - ^ 


/^(l) ( V''\ /^(l) ( V/>\ 

1 


- dss^G?^{X':)\sgniT-r), 


(69a) 

(69b) 

(69c) 


(70) 


where J dQ means to integrate over solid angle with unit 3-vectors fl, the 4-vector hi = 
(0, hi), the subscript R means the radial direction, and we dehne X" = x + (1/2) |T — T'\fl, 
X'J = x + {s/2)\T-r\Q, and 


|T-T'|/2 


g“ (.Y") = Gab(X") - Gab(x) = / drGAB.i(x + . 


( 71 ) 


The order-1 remainder term is 


Ri{T,T') = 


327r2 


dn{- 


1 r 


G?^(X") - Ggi(X") 

1 


- / *s"G®(A'")bgn(r-r'), 


(72) 

f3) 

where is the remainder after subtracting the second-order Taylor series. We can write 


GUx") = I 


^ rlT-T'l/2 


drGABjjxix + rQ) 


f 


T-T' 


V 2 




(73) 


When we apply the r and u derivatives from Eq. fl46|) . we can take {T, Z,X,Y) = 
{t,z,x,y) and calculate OIHq, O^Hq, d^Hi, d^Ro, and d^Ri. Applying u derivatives to 
Ho gives 


= 


dSvr^ 


Rtt,uu{x^ ^—) 


(74) 


For the derivatives with respect to r we have 


dlHG = 




R(x ), 


(76) 


and 


dlH^ = 


3207r2 


i^Rtt,tt{^) - - - aRii{x) + -R,tt{x) (3 -F In (-r_)) 


(76) 
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in the r —)• 0 limit. 

Applying u derivatives to i?o gives 

dlRo =^2] I 2 - GrrA^'")] 

dss'^Gtt,iAs')'^Q'sgaT, (77) 

where x'" = x + rQ and x"' = x + srQ. 

Now we have to take the second derivative of Ri with respect to r, which is T —T' in this 
case. This appears in three places: the argument of sgn in Eq. fl72ll , the limit of integration 
in Eq. ([73]), and the term in parentheses in Eq. flTHll . When we differentiate the sgn, we 
get 6{t) and 5'(r). but since ~ r^, there are enough powers of r to cancel the 6 or 
S', so this gives no contribution. When we differentiate the limit of integration, the term in 
parentheses in Eq. fl73]l vanishes immediately. The one remaining possibility gives 

d^Ri = j dQj dr I- [Gu,ijkA'") - Grr,ijkA'")] 

ds sgn r . (78) 




C. Derivatives with respect to ( 


To differentiate with respect to (, we must consider the possibility that x and x' are not 
purely temporally separated. We will suppose that the separation is only in the t and z 
directions and construct new coordinates (T, Z) using a Lorentz transformation that leaves 
X unchanged and maps the interval (T — T', 0) in the new coordinates to {t,() in the old 
coordinates. Then 

T-T'= sgur^r^^^, (79) 

and the transformation from (T, Z) to {t, z) is given by 


A = 


r C 

sgnr^— ^2 \ C x 


with the X and y coordinates unchanged. Then 


T-r 

0 


(80) 


(81) 


Now let M be some tensor appearing in i7(i). The components in the new coordinate 
system are given in terms of those in the old by 

MABC...=RMR'c---Mabc... (82) 


We would like to differentiate such an object with respect to ( and then set C = 0- The only 
place ( can appear is in the Lorentz transformation matrix, where we see 




= T 


.-1 


C=o 


0 1 
1 0 


13 


(83) 








and similarly, 


dlK\ 


= T 


.-2 


C =0 


1 0 
0 1 


(84) 


P = 


To simplify notation, we will define P and Q to be the matrices on the right hand sides. 
Reinstating x and y, 

/O 1 0 0\ 

10 0 0 
0 0 0 0 
\o 0 0 oy 

/ 1 0 0 o\ 

0 10 0 
0 0 0 0 
y 0 0 0 oy 

Now we can write the derivative of Mabc... as 


<2 = 


( 86 ) 


( 86 ) 


d^MABC... 




C=o 


(87) 


C=o 




Mabc... 
1 


C =0 

abc... 


= —[nVc ■■■ + ■■■ + ■■ OM.fc... = —pSa..Mafc... 

r_ 8 -^ r_ 


where p‘abc a rank-n matrix of O’s and I’s. With two derivatives, we have 


O'^Mabc... 


= dl{k\k^sk'c...)Mabc... 


C=o 


( 88 ) 


C=o 

Os ( 02 A ft \ xc 


+2(acA^)(acAi)5^ • • • + 2(acA^)5|(acA^) • • • + 
1 


M 


abc... 


C=o 


= -(Q»c--- + Wn^o--- + --- 

T_ 8 - 


xc 

"V" 

n 

ja^b Tjc 


+ P%Pi5^C--- + P15BPc--- + ---)Mabc.. 
1 


V 

{n—l)n 


— 2 ^ABC...^abc... 

T_ 

where ^ rank-n matrix of nonnegative integers. 

There are also places where T — T' appears explicitly in Hi. We can differentiate it nsing 
Eq. (179]), 


a^(T - T) 
dl{T-T') 


= 0 , 


C=o 


= —r 


.-1 


(89a) 

(89b) 


C=o 


14 

























Now we apply the operators and drdc^ to Hq, Hi, and -Ri. First we apply one C 
derivativ^ to Eq. fl69bp nsing Eq. fl87|l . 


dr 


C=o> 


' ^2 PttRab{x 


and two C derivatives nsing Eqs. (j 88 ll and (ISQah . 

1 


dlHo 


C=o 


dSvr^r^ 




(90) 


(91) 


Then we apply one C derivative to Hi, 

1 


dr I dc_Hi 


C=o> 


19207r2 


Pmt Rab,cd{x) - Pii ORab{x) + Rab{x) (In (-r_) + 2) 


and two ( derivatives to Hi, 
1 


(92) 




2=0 


6407r2 


pmfRab,cd{x) - ‘^RttA^) + OR{x) -^{q^OR^hix) 
rIquRA^) ) In (-r!) + ^ ( aRii{x) + ^R^x) 1 (1 + In (-r!)) 


(93) 

Finally we have to apply the derivatives to the remainder Ri. We can apply the 

(3) 

derivatives in two places, the Lorentz transformations and G\q. Since the three terms are 
very similar we will apply the derivatives to one of them 


Rc / cLVlgPAx + Y) 


C=o 


= / <in I + Y)d^Y‘ 


T 


(94) 


C=0y 


where we defined = (1/2)|T —T'|A“r2^. Then using Eqs. (|87)l and (I89al) . we find that that 
= (l/2)p“n* sgn r and taking into account the properties of Taylor expansions, 


d ^(3) 


dY^ 


G^Ax + Y) 


= G^ix"), 


(95) 


C=o 


(' 2 '] 

where G\^^ is the remainder of the Taylor expansion of Gab,c after subtracting the hrst-order 
Taylor series. 

Thus Eq. flUT]) becomes 


d^ I dnGfAx + Y) 


C=o 


da Rpf,G^2(x") + \G‘'Z(x'')p';a‘ 


sgn r 


(96) 


^ The Lorentz transformation technique we use here is not quite sufficient to determine the singularity 
structure of the distribution d^HQ at coincidence. Instead we can use Eq. (47) of Ref. to compute the 
non-logarithmic term in Hq for arbitrary x and x', which is then —Rab{x){x — x')°‘{x — a;')^/(487r2cr+). 
Differentiating this term gives Eq. (1901) and explains the presence of r_ instead of t in the denominator. 
The first term of Eq. arises similarly. 
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Using from Eq. fl75]) and 



Eq. fl96|) becomes 




dncPTix + Y)\<^=o = 



drGab,ijc{x + rn) 






(97) 




(98) 


-PiGtt,ija{x"') sgnr 




We could simplify further by using the explicit values of the p matrices, but our strategy here 
is to show that all terms are bounded by some constants without computing the constants 
explicitly, since the actual constant values will not matter to the proof. 

Applying the r derivative gives 


dr 



dVlG^^ix + Y) 




PttGab,ijk{x'") 


+ ^PiGtt,ija{x'") 




(99) 


We do not have to differentiate sgnr here, because the rest of the term is O(r^) and so a 
term involving S{t) would not contribute. 

The same procedure can be applied to all three terms. Terms involving X'J will get an 
extra power of s each time G is differentiated. The hnal result is 


dr I dcRiiT,r) 
1 


( 100 ) 


C= 0 y 

\r\/2 


327r2 


dVt 


dr 


'0 




^{pt-pTr)Gab,M")- / dssYuGab.Ms) 


~^{Gtt,ija{x"') — Grr,ija{x"')) — / dsS^p^Gu,ija{x'") 


sgn r. 


For two ( derivatives we can apply both on the Lorentz transforms, both on the Einstein 
tensor or one on each. 


aj / dQG^^ix + Y) 


/ 


C=o 


= / c^fi(^G-(x") + 


d^ 


G^\x + Y)d^Y‘^d^Y^ 


d ^(3) 


+2 


dY 
Ptt d 


G^^\x + Y)dlY 


2\ra 


C=0 


( 101 ) 


C-0 


GS>(j + y)acy' 


r dY’^ 


C=o 


Using Eqs. fl89bD and 




C=o 


2r 2r 2r * ’ 


( 102 ) 
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with h{ = qi - dj, while 




= 0 , 


(103) 


C=o 


since = 0 and hi* = 0. Using properties of the Taylor series as before, we can write 


32 




(b 


(104) 


so Eq. fllOip becomes 


[ dnGP^ix + Y) 

= [ dn 

J 

C=o J 


■ ^ab 1 

%gS>(x") + -ptpp^^;jx")n‘{V 


2 r 


2ptp‘GZ(x") + G‘S^(x")hl]n‘ 


(106) 


Using as in Eq. (ITT]) and and from Eqs. fl^ and fl73|) this becomes 


aj / dnGP^ix + Y) 


r /•hl/2 r /I ^ \ 2 

= J dfl J y2~ ]r[j ^ttGab,ijk{x'‘ 

Y-p-‘p)GuxM") + i\ ’’ 


J 


Yh^Gtt,ijk{x'") 


4 2\t\ 


2ptMGabM^'") 


(106) 


For all three terms 


d^RiiTX) 


C=o 


^,fdn f' " drU\-X 


2 r 


(107) 


(9,7 - qfAGatmix'") - / dss‘qtG^,j,(x":) 




4 2|r| 


:{Gtt,kab{x''') — Grr,kab{.x''')) — / dsS^G tt,kab{x'”) 


pMt - Pfr)Gab,M^"') + 7 ;K{GuM^'") - Grr,ljk{x'")) 


dss%2pM^G.b,kcK) + KGu,iA<')) 


sgnr. 


V. THE FOURIER TRANSFORM 


Eqs. dp, dZl, dZSl), (Hi, (Hi, (Hi, O, (EID, dli, dH, (Hni and dMl) include all 
the terms. To perform the Fourier transform we expand according to 
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Eqs. fH5|) and fH 6 |) and separate the terms by their r dependence, 


dudu'H-i + ^{ dlHo + ^iRo 


1 

2 ’ 


■{d^ + + 2drd(^) \Hq + Hi + -iRi 


T_ \TT^ / T_ T_ 

. 2 ' 


+ In (-r_)c2(t) + C3(t) + C4(t, r) 


( 108 ) 


where ci, C 2 , and C 3 are smooth and have no r dependence and C 4 is odd, Ci and bounded. As 
mentioned in Sec. IIV Al the functions yt depend on r but are smooth. Explicit expressions 
for the Ci are given in Appendix 

We now put the terms of Eq. f[108l) into Eq. (HOj) . and Fourier transform them, following 
the procedure of Sec. IV of Ref. |9|, to obtain the bound B in the form 


6 

B = 6-^^Bi. (109) 

i=0 

The hrst term in Eq. (llOSh is l/( 7 r^rl), and we proceed exactly as Ref. [^, except for the 
different numerical coefficient, to obtain 


Ro = 


1 

247r2 



( 110 ) 


Putting onlv Eq. fITT0|) into Eq. flT09|) rives the result for flat space. Fewster and Eveson na 
found a result of the same form, but they considered Tu instead of so the multiplying 
constant is different. Fewster and Roman [l^ found the result for null projection. Where we 
have 1/24, they had (n-£)^/12, where v is the unit tangent vector to the path of integration. 
Here v ■ i = R = l/((5\/2), from Eq. (02]), so the results agree. 

The remaining term requires more attention, because of the r dependence in yi. We 
write 

B,= r rfrGi(r) 4 e-*«" , (HI) 

Jo tt H 

with 

Gi{T) = j di{yi(t,T))g g (t+^'^ . (112) 

Then 

= 4g;"'(o) . (113) 

Applying the r derivatives to Gi gives 


GTir) 


dt 


r=0 -00 


d^ , 


r=0 


+ 3^2/1 (^,'r) 


12 /" 

r=0 

n(T\2 


- g'itr 


■g2/i(^(/"(^^(^ - 4:g'''{i)g{t) + 3g"{t) 


(114) 
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where the terms with an odd number of derivatives of the product of the sampling functions 
vanish after taking r = 0. 

Now Hi depends on r and i only through t + (A — l/2)r, so using Eq. fl6^ . we can 
write 

= ^ ^ ^ f d\{X-ll2)Yi(t”). (116) 

dr dr Jq dt Jq 

Then we integrate by parts and put all the derivatives on the sampling functions g, 


1 /.oo _ 


'0 


dX j^A - -j + Ag'{i)g'"{i) + g{t)g''"{t)) 

"1 ^ 


+3 dX I^A - Fi(f)(c/""(t)^(f) - g"{r> ) 


(116) 


Since we set r = 0, Ti has no A dependence and we can perform the integral. The result is 
^1 = y dtYi{rj{g''{tf - 2g'"{t)g'{f) + 2g''"{t)g{f)) . (117) 

For the term proportional to rT^, we have 


oo roo 


B, = 


71 


drG'2(r)^e-*«T 


(118) 


where 


G 2 {r) = j dty 2 {t,T)g g (t+^ 


We calculate this Fourier transform in Appendix [B] and the result is 

1 


Applying the derivatives to G 2 gives 

poo 

= dt 


G'"(r) 


r=0 J -00 


d? 


B; = gG;'( 0 ) 


3 d 


r=0 


9 {^ +2^2/2(t,r) 


{g'\t)g{^ - g\tr 


T=0 


(119) 


( 120 ) 


( 121 ) 


Again the only dependence of 1/2 on r is in the form of t” so we can integrate by parts 


POO /»! 


B 2 = —- / dt I dX 


' —00 J 0 


2 (^A - -J Y,it)i3g'ii)g"it) + git)g"\t)) 
+\ - y) ^ 2 (^(/'(^c/(^ - g''{i)g'{i)) 


and perform the A integrals 


-02 = i f dtY 2 (t)(g'(t)g"(t) - 3g"'(t)g(rj). 


( 122 ) 


(123) 
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For the term proportional to r_^, we have 


r°° (if 1 

Bs= - drG3(r)^e-‘«^ 


TT 


where 


G^(t) = 


cii(ci(t)+ 93(t,T))g (i-g (i +. 


(124) 


(126) 


Ref. calculated this Fourier transform, but the Fourier transform of 1/r^ given by Ref. 13 
was cited with the wrong sign in Eq. (105) of Ref. [^. The correct result is 




(126) 


Applying the derivatives to G 3 gives 

(f 


G'i{T) 


dt 


T=0 J -00 




2/3 (^,'r) 


d(^ + o(ci(^ + 2/3(^))(/(^)d(^ - d'(^ ) 


r=0 


(127) 


As before, we integrate by parts 


nOQ 1*1 

B3 = - di dX 

f J -00 Jo 


Integrating in A gives 

poo 




+ (1 - X)XYs{r)){9''{r)9{t) - c/'(t)2) 


(128) 


B. = 


dt 


c,{t)ig"it}git) - g'm + Y,it}{39"mi^ - 

15 


(129) 


The three remaining terms have Fourier transforms given in Ref. j^, so we hncH 

/ oo /* 00 

di / drg (i +—^ g (^i ——^lii\T\c 2 {i)sgiiT (130a) 

-oo J-oo ^ 2/ \ 2/ 

poo 

Rs = / dig{i)‘^{c3{i) + C5{t)) (130b) 


Rfi = 


dt 


TT 


— OO ^ —oo 


dr^g i+^)9(t- c^it, r), 


where we added 


C5(^ = -i2a + b)Ruu{i), 


(130c) 


( 131 ) 


which is the local curvature term from Eq. 0401) . 

The bound is now given by Eqs. ffTOOD . moj) . ffTTfl) . OT^ . OT^ . ffT30|) . 


Equation (I130a|) corrects an error of a factor of 2 between Eqs. (114) and (116) of Ref. 
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VI. THE INEQUALITY 


We would like to bound the correction terms Bi through Bq using bounds on the curvature 
and its derivatives. Using Eq. ([3]) in Eq. fl66ap . we hnd 


ym < 


3 


^max • 


(132) 


We can use Eq. (11321) in Eq. (I117p to get a bound on |i?i|. But will not be interested in 
specihc numerical factors, only the form of the quantities that appear in our bounds. So we 
will write 


Bl\ < jf\g]\x^\^Rra.. 


(133) 


where [g\ is an integral of some combination of the sampling function and its derivatives 
appearing in Eq. flll7p . We will need many similar functionals Jn\g\i which are listed in 
Appendix The number in the parenthesis shows the dimension of the integral, 



Similar analyses apply to B 2 and B^ and the results are 

1^21 < jf\9]\x^\Rm.. 
I-B3I < J3 \g\Rm 3 .-^ ■ 


(134) 


(135a) 

(135b) 


Among the rest of the terms in B there are some components of the form Rated,u which 
diverge after boosting to the null geodesic, as shown in Ref. j^. However we can show that 
these derivatives are not a problem since we can integrate them by parts. Suppose we have 
a term of the form 

roo poo 

Bn= di drLn{T,i) Rated, u{x), (136) 

J 00 J —00 

where Ln{T,i) is a function that contains the sampling function g and its derivatives. The 
u derivative on the Riemann tensor can be written 


Rabcd,u Rabed,t Ri 


^abed.v 


(137) 


The term can be reorganized the following way by grouping the terms with t and n, x, y 
derivatives 


poo poo 

Bn= di / drLnir, i){Af'^'^Rabed,tix) + A“R,fe,d,«(5;)), (138) 

J 00 J —00 

where are arrays with constant components and the subscript n denotes the term 

they come from. Here the greek indices a, fd, ■ ■ ■ = v,x, y. The term with one derivative on 
a can be bounded, while the term with one derivative on t can be integrated by parts, 

poo poo 

Bn = - di driL'^il T)Ai^^<^Rabed{x) + Ln{l T)A^^^^^Rabed,aix)) • (139) 

J oo J —oo 
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where the primes denote derivatives with respect to i. The sampling function is so 
L'{T,i) is still smooth and the boundary terms vanish. Now it is possible to bound this 
term, 


POO /*oo 


\Bn\< I dt I dr(|L^(t,r)|a®i?max + , 

where we dehned 


' OO «/ —OO 


qM = 


E 

abed a/3... 


A: 


abed a^... 


(140) 


(141) 


The same method can be applied with more than one u derivative. 

Now we apply this method to the integrals S 4 , and of Eq. fll30p . We start with 
Bi, which has the form 


S 4 = / dr In |r| sgn r / r) ( Rahcd,tt{x) + Rahcd,atix) 

, (142) 


where 


Lt{t,T) = g(t + T/2)g'(t - r/2) 


(143) 


After integration by parts 


Bi= I drill |r| sgnr / dt[ Ll{t,T)Af'^'^Rabcd{x) - L'^{t,T)Af'^'^°‘Rated,a{^) 


1 abeda : 


Taking the bound gives 


+L4(t,r)A“^i? 

abed^ap 


IB4I < E 

m=0 


(144) 


(145) 


Reorganizing B^ based on the number of t derivatives gives 


B. = 


' —OO 

poo 


dtL,{t} I Af<^RabcdA^) + Af^^RabcdM^^) + Af^<^^RabedM^) ) (146) 


= j dti^L'/i^Af^^RaUx) - L'^{t)Af^<^RabedA^) + L,{t)At^'^^Rated,a0{^) ) , 

where 

and the bound is 


^5(t) = d(^^ 


(147) 

(148) 


m=0 
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Finally the remainder term is 

/oo /*oo p pi ^ 

dt / drleilr) dn dx\ Af'^<^{X,n)RabcdM>^^) 

-oo J —OO J Jo L 

+A“(A, Q)R,bcciM>^^) + Q)RabcdM>^^) 

RAf^d^^^lX, hl)i?a6cd,a/37(Ahl)| sgnr , (149) 


where we changed variables to A = r/r and now arrays have components that depend 

on A and Al, and 

= 5'(^-'7'/2)fi'(t+ r/2). (150) 

After integration by parts 

/ CXD poo p pi ^ 

dr di dn dx\L'^iT,^At\x,n)RabM^) 

-oo J —oo J Jo L 

+L,{t, t)"A“(A, Q)RabcdA>^n) + L,{t, t)'A“^(A, Q)Rab,apAn) 
+Le{T,i)AA'^^'^{X,Q)Rab,a 0 -,{XXl) \ sgnr. (151) 


We dehne constants 


(m) 

«6 = 


E 

abed 0/3... 


da / dXAf^‘‘“'’"'(x,a) 


(152) 


and now we can take the bound 




m=0 


Putting everything together gives 


6 3 


B < +E E 4 ‘-”’is]r 


(m) 

max 


n=l 


n=4 m=0 


(153) 


(154) 


We can change the argument of the sampling function, writing g{t) = /(t/to)i where / is 
dehned in Sec. m and normalized according to Eq. CB, SO Eq. (I154p becomes 


dtTUn>m(tr > -^1 Ah f° */"(*/(»)" + ^ 


''0 ^OJ-to 

6 3 


n=l 


EEj<‘-”)|/]flW(; 


m+2 

0 


n=4 m=0 


(156) 


where we used jA[g] = tQ^Jn\f]. We can simplify the inequality by dehning 

^ = / f'iafda = f'\t/tofdt , 


(156) 
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( 157 ) 


and 


n=4 

pin) ^ J^4-n)jy| _ 


Then Eq. (I155p becomes 

j dtT^u{w{t))g{tf > 

5-2 I 1 


F+Y. + I] 


to 1 247r2 

^ ^ m =0 

We will nse this resnlt to prove the achronal ANEC. 


n=l 


(158) 


(159) 


VII. THE PROOF OF THE THEOREM 


We nse Eq. fll59p with w{t) = ^viv^t) and integrate in r] to get 

/ m po 

dr] Tuui^viv,t))f{t/tof> (160) 

Vo J—to 

_j]p_ I ^ F I p(m) p(m) y.m+2 , |-S|3-np(n)p +>^-11 

^2^3 \ 2471-2 -“-max^^O + P I ^ -O-max^o ^ • 

0 1^ m=0 n=l ) 

As 5 —)■ cx), to —t 0 bnt F^'^\ F^'^\ i?max, and R ^X are constant. Now x“ = x“/5, and using 
Eqs. m, (USD, m , |x“| < Ui + \/25to- Thus as 5 —>■ oo, x“ ^ 0. Therefore only the hrst 
term in braces in Eq. fll60p survives, so the bound goes to zero as 


Vo 

6Hl 


j:2a—l 
~ 0 


(161) 


Equation fll60p is a lower bound. It says that its left-hand side can be no more negative 
than the bound, which declines as But Eq. flT^ gives an upper bound on the same 

quantity, saying that it must be more negative than —Ato/2, which goes to zero as to ~ 5“". 
Since a < 1/3, the lower bound goes to zero more rapidly, and therefore for sufficiently large 
5, the lower bound will be closer to zero than the upper bound, and the two inequalities 
cannot be satished at the same time. This contradiction proves Theorem 1. 

The ambiguous local curvature terms do not contribute in the limit rjo ^ oo because they 
are total derivatives proportional to 



drjR^uuix)) = 0 . 


(162) 


VIII. CONCLUSIONS 

This work completes the proof of ANEC in curved spacetime for a minimally coupled, 
free scalar held, on achronal geodesics traveling through a spacetime that obeys NEC. The 


24 














techniques are similar to those of Ref. [3|, but that paper required an unproven conjecture. 
Here, we use the general absolute quantum inquality of Fewster and Smith to derive a 
null projected quantum inequality, slightly different from our previous conjecture, and use 
that inequality, Eq. fllSSp . to prove achronal ANEC. Equation 01551) h as the form of the flat- 
space null-projected quantum inequality of Fewster and Roman [l2(], plus correction terms 
which vanish as one considers more and more highly boosted timelike paths with smaller 
and smaller total proper time in the limiting process above. 

The result of this paper concerns integrals of the stress-energy tensor of a quantum held 
in a background spacetime; we have so far not been concerned about the back-reaction of 
the stress-energy tensor on the spacetime curvature. This analysis is correct in the case 
where the quantum held under consideration produces only a small perturbation of the 
spacetime. Thus we have shown that no spacetime that obeys NEC can be perturbed by 
a minimally-coupled quantum scalar held into one which violates achronal ANEC. Thus no 
such perturbation of a classical ^acetime would allow wormholes, superluminal travel, or 
construction of time machine^ [l|. 

What possibilities remain for the generation of such exotic phenomena? One is that the 
quantum held is not just a perturbation but generates enough NEC violation to permit 
itself to violate ANEC also. We argued against this idea on dimensional grounds in Ref. j^. 
Another is that there is a held that violates NEC but obeys ANEC, and a second held, 
propagating in the background generated by the hrst, that violates ANEC. This three-step 
process seems unlikely but is open to future investigation. 

There is also the possibility of diherent helds. We have not studied higher-spin helds, but 
these typically obey the same energy conditions as minimally-coupled scalars. Of more inter¬ 
est is the possibility of a non-minimally coupled scalar held. Such helds can produce ANEC 
violations even classically [l^, T3| with large enough (Planck-scale) held values. However 
these situations seem unphysical since the ehective Newton’s constant becomes negative as 
the held value increases. In the case of a wormhole [l^, the ehective Newton’s constant 
must be negative not only inside the wormhole but in one of the asymptotic regions. If one 
disallows Planck-scale held values, there are restrictions on non-minimally coupled classical 
(l^ and quantum 201 helds, but these restrictions are not in the form of the usual quantum 
inequalities. Whether there is a self-consistent achronal ANEC for non-minimally coupled 
scalar helds remains an open question. 
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^ More specifically, ANEC rules out compactly generated causality violation [l^. A causality violating 
region is compactly generated if the generators of the Cauchy horizon followed into the past enter and 
remain within some compact set. But Ori 1^ argues that one can ensure that closed causal curves develop 
even without compact generation, by arranging the formation of a null hypersurface spanned by closed 
null geodesics at the future boundary of the domain of dependence. 
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Appendix A: Explicit expressions 


Here are the explicit expressions of the functions q used in Eq. (11081) : 

Cl = ^ ( - R{X) + [pt - 

C2 = ( - ^RUx) + {vt +1^) ^Ra,{x) + ^ RA^)^ 

C3 = ggQ^2 “ 2 Rab,cdA + + 2 ^,tA) 

+PiARabix) + ^PuR,ab{x)^ 


'hl/2 


C4 = 


2567r2 


dr j dQ^dl!^^[Gu 4 x'") - Grr,iA')] - ds s^Gu,i{x"s') 


i [Gtt.jkA'') - Grr.jkA")] f dss^Gu,ijkix':') 


^(Pu - Pfr)Gab,,kA") - / dsA^Aab,Ax'n 


I 1 

+ (1-7^ 

+ f(Gu,jka(x"') - Grr,Jka(x"')) " / ^ 


, 1 r 
+2 --^ 


fr 


2 (du-C)Gab,ykA")- dssAtGabMA) 


-.Pip'j 


■,{Gtt,kabix''') ~ Grr,kab{x'")) 


dss^GttMbix's) 


r 


\T 

^ 3 /r)„c„ab 


pMApfr)Gab,kcA'') + 2h\iGu,l,kA") - Grr,ljkA")) 


- / dss\2p'rp'^,AabMA) + KGtt,ijki<))\^'^ 


And here are the integrals of the sampling function; 


/ OO 

c?^(aii|^""(^)|^(^) + ai2\9"'{t)g'{t) \ + ai3/(f)^) 

-00 

/ OO 

dt(a2i|/'(t)|^(t) + a22\9"{t)g{t)\) 

■OO 

/ OO 

c?^(a3i|^"(^)|^(^) + a32^'(^)^) 

-OO 

/ OO poo 

dt / dt' |ln \t - t'W iaii\g'''{t')\g{t) + a42|/(t)5''(f')l) 

-OO J —OO 

/ OO poo 

dt / dt'|ln |t - f'll (a43|fi'"(t')lfi'(^) + «44|5''(f)5''(f')l) 

-OO J —00 

/ OO poo 

dt / dt'|ln |t - t'll a45|5('(t')|5((t) 

-OO J —OO 


(Ala) 

(Alb) 

(Ale) 


(Aid) 

(A2a) 

(A2b) 

(A2c) 

(A2d) 

(A2e) 

(A2f) 
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4“’Is] = 

poo 

/ dt{a^i\g''{t)\g{t) + a^ 2 g'{tf) 

/ —OO 

(A2g) 

4”’lsl = 

poo 

/ dt a^^\g'{t)\g{t) 

J —OO 

(A2h) 

Jt’ls] = 

poo 

/ dta^ 4 ,g{tf 

J —OO 

(A2i) 

4%] = 

poo poo 

/ dt / dt'{a,,\g'"itMt') + a, 2 \ 9 "it)g'it')\) 

J —OO J —00 

(A2j) 

4“ la] = 

poo poo 

/ dt dt'{ae 3 \g"{t)\g{t') + aQi\g'{t)g'{t')\) 

J —OO J —OO 

(A2k) 

4'‘’l9l = 

poo poo 

/ dt dt' aG 5 \g'{t)\g{t') 

J —OO J —00 

(A21) 

1 

poo poo 

/ dt dt' aGeg{t)g{t'), 

J —OO J —00 

(A2m) 


where ank are positive constants that may depend on a^\ 


Appendix B: Fourier transform 

We follow the procedure of Ref. to calculate 

B2= ( —[ (Bl) 

Jo J —oo 

where 

G 2 {t) = j dty 2 {t,T)g g . (B2) 

and 

= (B3) 

This is the Fourier transform of a product so we can write it as a convolution. The function 
S 2 is real and odd, so its Fourier transform is imaginary, but G 2 is also real and odd, thus 
the Fourier transform of their product is real. We have 


Bo = 


27r2 


/ dCG^i-^ - OUC) • 


(B4) 


We can change the order of integrals and change variables to g = — ( which gives 

1 


Bo = 


27r2 

1 


/ dC dgG2{g)h{,C) 

' —00 J c, 

/ oo PT] 

/ dCkiO- 


(B5) 


The Fourier transform of S 2 is 13 


hiC) = -inC^QiC) , 


(B6) 
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and 


•V 


ITT 


From Eq. flBSp we have 


52 = -—/ dr]G2{vW. 

OTT ./n 


Using /'(O = we get 


(B7) 

(B8) 


^2 


I r°° ^ 

dvG'^v)- 


(B9) 


The fnnction G 2 is odd bnt with three derivatives it becomes even, so we can extend the 
intergal 


B 2 


1 



dr]G'”(7]) 


1 

6 


G'"(0). 
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